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Abstract 

We study classical chaotic motions in the Berenstein-Maldacena-Nastase (BMN) 
matrix model. For this purpose, it is convenient to focus upon a reduced system 
composed of two-coupled anharmonic oscillators by supposing an ansatz. We examine 
three ansatze: 1) two pulsating fuzzy spheres, 2) a single Coulomb-type potential, and 
3) integrable fuzzy spheres. For the first two cases, we show the existence of chaos by 
computing Poincare sections and a Lyapunov spectrum. The third case leads to an 
integrable system. As a result, the BMN matrix model is not integrable in the sense 
of Liouville, though there may be some integrable subsectors. 
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1 Introduction 

One encounters non-linear differential equations everywhere in theoretical physics. String 
theory is no exception, as a matter of course. For example, when a string is moving on a 
curved space, the classical equation of motion becomes non-linear generally. If the curved 
background has a special coset-structure like AdSs x S^, the string theory is classically in- 
tegrable in the sense that the Lax pair exists [^. Then classical solutions can be studied 
analytically, for example, by using the classical inverse scattering method. However, in 
most cases, one cannot solve analytically non-linear differential equations. The systems are 
non-integrable and exhibit chaos. 

Recently, chaotic string solutions have been well studied in various ways |2]-[^ . Motivation 
lies on potential applications to the AdS/CFT correspondence (^-[^ beyond integrability [lo) , 
but it has not been clarified yet what corresponds to chaotic strings in the gauge-theory 
side. Note that chaotic motions cannot be described with analytic functions. This is a 
characteristic of chaos. So it may be intriguing to ask whether the correspondence holds 
even in the non-analytic region or not. One can test it at a deeper level. 
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Along this direction, one may focus upon the M-theory dynamics. A matrix model 


description of light-front M-theory is proposed by Banks-Fischler-Shenker-Susskind 11 . It 


is called the BFSS matrix model. Chaos in this model was studied in [1^ by following 
the way in classical Yang-Mills (YM) theories (^. However, there are some unsatisfactory 
points. In the BFSS matrix model, there are flat directions which may cause unbounded 
motions. Then the motions may conflict with the boundedness of chaos during long-time 
numerical computation of classical trajectories. Another point is a technical problem. In (I^ , 
the largest Lyapunov exponents were computed by reducing the system to simple models. 
However, it does not seem that the results converge on deflnite values. In addition, the 
time scale of (would-be) convergence is also much larger than the inverse of the Lyapunov 
exponent. Thus, it would be worth revisiting chaotic motions in the BFSS matrix model. 


A more appropriate system is the Berenstein-Maldacena-Nastase (BMN) matrix model 
(l4| . It contains mass term^ and flat directions are lifted up. Hence there is no subtlety 
with the boundedness of chaos. In addition, the Myers term may give rise to a separatrix 
structure in the phase space. It means a bifurcation of the behavior of classical solutions and 
may be a source of chaos. In total, the BMN matrix model is a good laboratory to study 
chaotic motions. 


In this paper, we will study classical chaos in the BMN matrix model. For this purpose, 
it is useful to reduce the full system to a simple model composed of two-coupled anharmonic 
oscillators by supposing an ansatz. We examine three cases: 1) two pulsating fuzzy spheres, 
2) a single Coulomb-type potential, and 3) integrable fuzzy spheres. For the first two cases, 
we show the existence of chaos by computing Poincare sections and a Lyapunov spectrum. 
The third case leads to an integrable system. As a result, the BMN matrix model is not 
integrable in the sense of Liouville, though there may be some integrable subsectors. 


This paper is organized as follows. In section 2 we give a brief review of chaos. In 
particular, the characteristics of chaos in the Hamiltonian system. Section 3 considers the 
Henon-Heiles system as an example of non-integrable Hamiltonian systems. In section 4 we 
consider chaotic motions and an integrable subsector in the BMN matrix model. Section 5 is 
devoted to conclusion and discussion. In Appendix A, the computation scheme of Lyapunov 
spectrum is summarized. 


^ For chaotic motions in a YM-Higgs system, see 


15 
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2 What is chaos? 


Probably, anybody has said the term, “chaos”. For example, when you see the top of the 
desk in a mess, you would say “It’s chaos”. So it seems likely that “chaos” implies that some 
things are completely messed up. In fact, this point captures an aspect of chaos, randomness. 

There would be no universally accepted definition of chaos. The following is a fist of 
characteristics of chaos commonly seen in almost any books: 

1. High sensitivity to initial conditions, 

2. Randomness of motions (continuous spectrum), 

3. Non-periodic and bounded motions (non-analyticity). 

The item 1) is significant and it is sometimes called the Butterfly effect. One may know 
the following famous sentence, “Does the flap of a butterfly’s wings in Brazil set off a tornado 
in Texas?” Here the flapping wing corresponds to a small change of the initial condition and 
it may cause a chain of events leading to large scale phenomena. 

Suppose a continuous dynamical system described by a set of non-linear differential equa¬ 
tions. The system is deterministic because the dynamics is determined by fixing an initial 
condition. But one cannot predict the motions generally after long time has passed because 
of the non-linearity. That is, it is impossible to forecast weather in the far future. This indi¬ 
cates that some information is lost dynamically and it is measured by the Kolmogorov-Sinai 
(KS) entropy. This entropy production is concerned with the randomness [the item 2)]. 

The item 3) can be intuitively understood. If the motion is periodic, then it does not 
satisfy the other properties. It exhibits the line spectrum and no sharp sensitivity to initial 
conditions. To see that the boundedness is necessary, let us consider a simple system x = x. 
The orbits run away exponentially to infinity, hence it seems to exhibit high sensitivity to 
initial conditions. However, it just means an attracting fixed-point at infinity and does not 
indicate chaotic motions. The recurrence of motions should be required implicitly so as to 
exclude fixed points. Non-analyticity follows from the continuous spectrum. The smooth 
behavior in the Fourier space indicates a random, non-analytic motion in the real space (i.e., 
the motion is continuous but not differentiable everywhere). Hence chaotic motions do not 
behave in an analytic way in comparison to (quasi-)periodic ones. Indeed, if the motion is 
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described by an analytic function, then one can predict the motion for an arbitrarily long 
time and it does not satisfy the item 1). 

There are three methods to see chaos, 1) Poincare sections, 2) Lyapunov spectrum and 
3) power spectrum. Since some readers are not familiar to these issues, it would be helpful 
to introduce each of them below. 

2.1 Poincare sections 

To study complex trajectories in the phase space, it is helpful to use a Poincare section, which 
is a surface of section intersecting the flow in the phase space. It is deflned so that all of the 
possible orbits cross the section repeatedly. The dimension of the flow (i.e., the dynamical 
system) is the number of the first-order ordinary differential equation^ and the dimension 
of Poincare sections is one lower than it. On a Poincare section, an orbit is mapped to dots, 
which are the intersection points of the orbit and the Poincare section. Since the dynamics 
is governed by a set of differential equations, the original trajectory can be reconstructed 
from the dots on the Poincare section. In this sense, the flow and the Poincare section are 
equivalent, while there is a degree of freedom to choose the surface in practice. It is the usual 
to pick up only the dots at which the trajectory intersects from only the one side to the other 
side. This constraint allows any dots on a Poincare section to uniquely determine the next 
point of the trajectory in autonomous Hamiltonian systems. Suppose a Hamiltonian system 
with a four-dimensional phase space (qi, g 2 ;Pi,P 2 ) • A Poincare section is a two-dimensional 
surface, for example, with ^2 = 0 and p 2 > 0. 

By changing initial values, a set of trajectories is generated and it contains several kinds 
of motions. The Poincare section is helpful to figure out their behaviors. Periodic orbits 
are projected onto the Poincare section as flnite numbers of dots on closed contours. A 
periodic trajectory comes back to a point on the Poincare section after its period. Quasi- 
periodic orbits are also projected as dots on closed contours, but the resulting portrait 
becomes dense. They are understood as non-resonant tori. Completely chaotic orbits result 

^ Note that, in an autonomous Hamiltonian system, all of the orbits are restricted onto a lower dimensional 
manifold due to the constancy of the energy. In this sense, the dimension of this dynamical system is one 
lower than twice the dimension of real space. In general, the dimension of a dynamical system is lower than 
twice the dimension of real space by the number of extra conserved charges. 
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in randomly scattered dots on the Poincare section. These behaviors depend on initial values 
and generally chaotic motions coexist with (quasi-)periodic motions on the section. 

Suppose that a small perturbation makes an integrable system be non-integrable. Then 
the Kolmogorov-Arnold-Moser (KAM) theory [I6| - [T8| claims that trajectories incline to be 
chaotic gradually as the perturbation becomes larger. When the perturbation is small 
enough, most of the orbits are (quasi-)periodic. These are called KAM tori, which are 
the remnants of Liouville tori in the unperturbed integrable system. In other words, almost 
all of the Liouville tori survive the small perturbation, though the tori are distorted and 
translated. As the perturbation becomes larger, chaotic trajectories appear and coexist with 
KAM tori. Then the region of chaotic trajectories enlarges while the KAM tori tend to get 
broken. When the perturbation is larger, KAM islands and islets are in the sea of chaos. 

2.2 Lyapunov spectrum 

The next is to introduce the Lyapunov spectrum. It is very useful to measure the dynamical 
information loss. 

Let us first introduce a Lyapunov exponent. We basically follow the description in [I^ . 
Suppose that a solution x(t) is chaotic and consider a nearby point x(t) + 5{t). Then the 
deviation 5{t) grows exponentially like 

||h(t)||^||h(0)||e^T (2.1) 

Here A is a positive constant called a Lyapunov exponent. This behavior indicates that the 
deviation of the initial value 5(0) is amplified exponentially even if it is very tiny. This is 
nothing but a realization of sensitive dependence on initial conditions. 

In general, the phase space is higher-dimensional. Then there are N different Lyapunov 
exponents for an A^-dimensional phase space. An intuitive explanation is the following. 
Consider time evolution of an infinitesimal sphere of perturbed initial configurations. During 
its evolution, the sphere tends to be distorted into an infinitesimal ellipsoid. Let Skit) [k = 
1,..., A^) denote the length of the fc-th principal axis of the ellipsoid. Then the deviations 
behave as 

Skit) ~ 5fc(0) , Afc : Lyapunov exponents. 
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The set of the exponents generates a spectrum called the Lyapunov spectrum. For large 
t , the diameter of the ellipsoid is controlled by the most positive , and it is called the 
largest Lyapunov exponent. 

The Lyapunov spectrum is also related to the dynamical entropy production. For ex¬ 
ample, the sum of all of the positive Lyapunov exponents is identified with the KS entropy 
when the system is ergodic. It would be very interesting to look for an application of the 
KS entropy in the context of string theory and M-theory. 


2.3 Power spectrum 


The power spectrum is defined as 


I(oj) = lim — 

T—)-oo 27 r 




( 2 . 2 ) 


where {y{t)y{s)) is a two time correlation function. It is quite helpful when one needs to 
distinguish quasi-periodic motions and chaotic motions. When the ratio of periods is irra¬ 
tional, the former leads to dense orbits in the phase space (the denseness). On the other 
hand, the latter is related to a stochastic motion (the randomness). In the Hamiltonian 
system, there is no attractor because the phase space volume has to be preserved due to 
Liouville’s theorem. Hence it is often difficult to distinguish them at a first glance. How¬ 
ever, the power spectrum shows clear distinction. For quasi-periodic motions, the spectrum 
contains sharp peaks associated with periods (the line spectrum). On the other hand, the 
continuous spectrum appears for chaotic motions. 


In the following, we will focus upon the first two methods, Poincare sections and a 
Lyapunov spectrum. 


3 Warming up — the Henon-Heiles system 


In this section, let us consider the Henon-Heiles system 20 . This is a famous example of 


non-integrable Hamiltonian systems containing chaotic motions. This model was originally 
proposed by Henon and Heiles (1964) in the study of celestial mechanics. 
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Henon-Heiles system 

The Hamiltonian is given by 


H = ]^{p\+pI) + ^(^1 + <il) + « (^^1^2 - , (3.1) 

where a is an arbitrary real constant parameter. The system is composed of two particles 
with the same mass (normalized to 1) and the interaction is third-order. The positions of 
the particles are described by qi and q 2 , and the conjugate momenta are pi and p 2 ■ It is 
well-known that this system is non-integrable. 


Before going to the detailed analysis, it may be interesting to see that the form of the 


interaction term in (|3.1|) is very important. If the interaction term is slightly modified like 

1 


a 1 qlq2 - -ql 


a ( qlq2 + -^2 


then the system becomes classically integrable. The resulting system is often called “anti 
Henon-Heiles system”. 


The system with (3.1) is equivalent to a periodic lattice with three particles, 


H — Hq + -a [(Ql — (52)^ + (Q2 ~ Qs)^ + (Qs ~ Qi)^] ) 

Hq = -^{Pi + P2 + P3) + 2 [(^1 “ ^2)^ + (Q2 — Qs)^ + (Qs — Qi)^] ■ (3.2) 


The unperturbed Hamiltonian Hq describes a system of three connected harmonic oscillators. 
Hence the equations of motion for Hq are exactly solved by using the normal coordinates. 
By using the normal coordinate for Hq , one can rewrite the total Hamiltonian H. Then, 
after dropping the total momentum, the Henon-Heiles system is reproduced. 


Here one may sum up higher-order interaction terms with appropriate numerical coeffi¬ 
cients and obtain an exponential-type potential. The resulting system is nothing but a Toda 
lattice system [^. Hence, inversely speaking, the Henon-Heiles system may be regarded as 
a truncation of the Toda lattice system. 


Poincare sections and Lyapunov spectrum 

Let us see chaotic motions in the Henon-Heiles system ilp. It is a good exercise to under¬ 
stand chaotic motions in the Hamiltonian system. Here we will consider Poincare sections 
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02 02 
(a) Poincare section with E = 0.0833 (b) Poincare section with E = 0.1250 



% 


(c) Poincare section with E = 0.16667 


Lyapunov spectum (E=0.16667) 



t(time) 

(d) Lyapunov spectrum with E = 0.16667 


Figure 1: Poincare sections and a Lyapunov spectrum. Different trajectories are indicated by dots with 
different color, depending on initial conditions. 


and a Lyapunov spectrum. These quantities can be evaluated numerically and the results 
are shown in Fig.[T} 

Figures (a), (b) and (c) are Poincare sections for E = 0.0833, 0.1250 and 0.16667, 
respectively. The sections are taken with gi = 0 and gi > 0. We set that o = 1. In Fig.[^ 
(a), most of the solutions are (quasi-)periodic (KAM tori). In Fig.[^ (b), chaotic motions 
appear but still (quasi-)periodic motions coexist. In Fig.[^ (c), most of the solutions are 
chaotic. 

Figure (d) is a Lyapunov spectrum with E = 0.16667. It is computed by taking 
an initial condition randomly from the region around qi = 0.0, pi = 0.35, q 2 = —0.115, 
P 2 = —0.443 with the radius 5.0 x 10“^. The largest Lyapunov exponent is Ai ~ 0.15. 
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This positive value indicates high sensitivity to initial conditions and provides a positive 
support for the existence of chaos in this system. As a feature of the Hamiltonian system, 
the spectrum is symmetric and the total sum of all of the exponents becomes zero. The 
remaining two exponents become zero due to the energy conservation. 

The features presented in Fig.are characteristics of chaos in the Hamiltonian system. 
In the next section, we will show similar results in the case of the BMN matrix model. 


4 Chaotic motions in the BMN matrix model 

We study here classical chaotic motions in the BMN matrix model. It is useful to reduce the 
full theory to a simple model so as to extract a mechanism which induces chaotic motions. 
We examine three reduction ansatze and argue whether chaotic motions exist or not. 


4.1 Setup 


The BMN matrix model 14 was originally proposed as a realization of M-theory on the 


maximally supersymmetric pp-wave background 22 . That is, it can be regarded as a de¬ 


formation of the BFSS matrix model fTTl. This model can also be derived as the matrix 


regularization of a supermembrane theory on the pp-wave [23, 24 by following a seminal 
work (^. 


The BMN matrix model is a one-dimensional U{N) gauge theory composed of matrix¬ 
valued variables X^{t) (r = 1,... ,9), a gauge field A and 16 fermions. For the present 
purpose, we will focus upon the bosonic sector. 

The bosonic part of the classical action is given by 


/ 


dtTr 


l(AVy + i[X'-,V|" 






(4.1) 


where i,j = 1, 2,3 and o, 6 = 4,..., 9 . The covariant derivative Dt is dehned as 


DtX^ = X^ -i[A,X^ 
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The symbol denotes the derivative with respect to t. The real constant parameter /i 
measures the deformation. When /i = 0, the system is reduced to the BFSS matrix model. 


When /r 7 ^ 0, mass terms and the Myers term 26 are turned orj^ 


In the BMN matrix model, the classical trajectories are bounded due to the mass terms. 
Hence there is no flat direction in comparison to the BFSS case. This characteristic would 
be so important as to study chaotic behavior definitely for sufficiently long time. 


Due to the presence of the Myers term, there exist fuzzy sphere vacua as well as the trivial 
vacuum X'^ = 0 . Thus the BMN matrix model has the rich vacuum structure. Furthermore, 
the Myers term leads to a separatrix structure in the phase space and may provide a new 
source of chaotic motions. 


We will work with the Weyl gauge A = 0 hereafter. Then the classical equations of 
motion are given by 

9 2 ^ 

0 = W + [X^ W]] + (I) + {i = 1,2,3), 


r=l 

9 


j,k=l 


0 = W“ + 5][X^[X^X“]]+(^) (a = 4,...,9). 


r=l 


(4.2) 


In addition, one has to take account of the Gauss law constraint, 

9 




(4.3) 


r=l 


This comes from the equation of motion of A . 


In the next subsection, we will study classical solutions of the equations of motion (4.2) 


equipped with the Gauss law constraint (4.3). 


4.2 Reduction to simple models 


To argue classical chaotic motions in the BMN matrix model, it is sufficient to study a simple 
reduced system by imposing an ansatz. We will examine three examples of the reduction 
ansatz below. 


^ Note that the mass ratio for X® and X“ is fixed by requiring the maximal 32 supersymmetries for the 
pp-wave background T4p2 . The coefficient of the Myers term is fixed as well. In this sense, supersymmetries 


are implicitly concerned with our analysis, though the present argument is restricted to the bosonic sector. 
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(1) Two pulsating fuzzy spheres {N = 2) 


First of all, let us suppose the following 2x2 ansatz of classical solutions: 

X* = rit) j, = x(t) ^ (a' = 4, 5, 6), 

= X^ = X^ = 0. 


(4.4) 


Here cr* (i = 1,2,3) are the standard Pauli matrices. Then x(t) and r(t) are real functions 


to be determined. This ansatz (4.4) describes two pulsating fuzzy spheres and satisfy the 


Gauss law constraint (4.3). 


With the ansatz (4.4), the classical equations of motion (4.2) are reduced to 

2 


0 = r + 


r — jdr^ + 2r^ + 2rx^ , 


fX 


0 = X + ( ^ ) X + 2x^ + 2xr^ . 


(4.5) 


The resulting system (|4.5|) is a non-linear dynamical system in which two anharmonic oscil- 

also be 
1 


lators are coupled. Note that the system (4.5) can also be obtained from the Lagrangian, 

1 


T ^-2 ^ ^ 2 ^4 22 

L = -r — -r r — — 4- -x — - — x — -x — r x 
2 2 V 3/ 2 2 V6/ 2 


(4.6) 


Poincare sections are presented in Figs.[^ (a)-(e). The sections are taken with r = 0 
and r > 0. For E = 0.001, 0.01, 0.1,1.0 and 10 with jx = 2.0, the sections are numerically 
computed. When E = 0.001, the energy is sufficiently small and only (quasi-)periodic 
trajectories appear [Fig.[^ (a)]. When E — 0.01, one can see chaotic motions between two 
layers [Fig.[^(b)]. This is a typical example of local chaos. When E = 0.1, KAM tori around 
the outer circle are destroyed and a ring of global chaos appears, according to the Chirikov 
overlap criterion. There is a big island in the sea of chaos [Fig.[^ (c)]. When E = 1.0, the 
island in Fig.[^ (c) splits into two islands [Fig.|^ (d)]. When E — 10, (quasi-)periodic orbits 
appear outside the sea of chaos [Fig.[^ (e)]. This is the onset of re-ordered motions. This 
phenomenon is quite natural because motions tend to be harmonic oscillation as the energy 
becomes much higher. 

The Lyapunov spectrum for FI = 10 and = 2.0 is presented in Fig.[^(f). It is computed 
by taking an initial condition randomly from the region around r = 0.0, Pr = 3.8, x = —0.50, 
P 2 = 2.33 with the radius 5.0 x 10“^. The largest Lyapunov exponent is approximately 0.302 , 
which is signihcantly non-zero. This result gives a support for the existence of chaos. 
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Poincare secton (E=0.001, n=2.0, at r=0.0) 


Poincare section (E=0.01, n=2.0, at r=0.0) 





(a) Poincare section with E = 0.001 


(b) Poincare section with E = 0.01 


Poincare section (E=0.1, ^=2.0, at r=0.0) 



Poincare section (E=1.0, ^=2.0, at r=0.0) 



(c) Poincare section with ^ = 0.1 


(d) Poincare section with ^ = 1.0 



(e) Poincare section with ^ = 10 


(f) Lyapunov spectrum with E = 10 


Figure 2: 


Poincare sections and a Lyapunov spectrum with the ansatz (4.4). 
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(2) A single Coulomb potential {N = 3) 

As the next example, let us consider the following 3x3 ansatz: 


/2/3 


= ^/6x{t) 


= ^/Qx{t) 


\ 


-1/3 


/-1/3 


/-1/3 


, = \/Qx{t) 


-1/3/ 

\ 


\ 


2/3 




-1/3/ 


-1/3 




2/3/ 




\ 



^ 0 

ze*‘^\ 

II 

0 —ze“ 


II 


0 



0 

) 





0 / 


( 0 -ze-*^ \ 




\ 


= yit) 

0 


, X^ = y{t) 

0 e“ 


5 



V 



e*^ 0 / 



i 0 6*“^^ 



( 0 

e-A \ 



= y{t) 

0 

1 

A'’ = y(t) 

e*^ 0 




0 / 




0^ 




(4.7) 


Here x{t) and y{t) are real functions to be determined. The angle variable ip{t) is cyclic and 
the associated constant of motion is L = y^ip. This ansatz (4.7) was originally proposed by 


Arnlind and Hoppe (^. Then the resulting differential equations are given by 


0 = x + X + 12xy‘^, 




0 = y+ (^0 y + 8y^ + 12x‘^y - — . 


The system (4.8) can also be obtained from the Lagrangian, 
1 


r --2 o f 2 . ^ -2 ^ f 2 q 4 ^ « 2 2 

2 \qJ 2^ 2 Ve/ ^ ^ 2u2 ^ 


2y2 


(4.8) 


(4.9) 


This system can be seen as a system composed of two-coupled anharmonic oscillators with 
a repulsive Coulomb-type potential. 

Poincare sections are shown in Figs.|^ (a)-(c). The sections are taken with a; = 0 and 
x > 0. For E = 2.0, 2.8 and 5.0 with fi = 3.0, the sections are numerically computed. 
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When E = 2.0, trajectories are (quasi-)periodic [Fig.[^ (a)]. The left KAM tori are pushed 
by the right KAM tori, but the energy is not suihciently high and hence there is no chaotic 
behavior. When E = 2.8, the energy is high enough to destroy the contiguous KAM tori in 
Fig.|(a). Hence chaotic motions appear along the contiguous line [Fig.[^(b)]. This behavior 
would be understood as local chaos again, according to the Chirikov overlap criterion. When 
E = 5.0, almost trajectories become chaotic and few ones remain (quasi-)periodic. 

The Lyapunov spectrum in Fig.|^ (d) also supports the existence of chaos. It is computed 
by taking an initial condition randomly from the region around x = 0.95, = 0.0, y = 0.5, 

P 2 = 1.44 with the radius 5.0 x 10“^ . 



PoincatB section (E=5.0, L=1.0, n=3.0 at x=0.0) 



Lyapunov spectrum (E=5.0,L=1.0, p =3.0) 



(c) Poincare section with ^ = 5.0 


(d) Lyapunov spectrum with E = b.Q 


Figure 3: Poincare sections and a Lyapunov spectrum with the ansatz (4.7). In Fig. (a), two KAM tori 


are contiguous. In Fig. (b), local chaos appears around the boundary. 
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(3) Integrable fuzzy spheres {N = 4) 

Let us concentrate on the SO{3) snbsector and consider the following ansatz: 


X* = - 


1 ( r{t) cr* x(t) cr* 
x{t) cr* r{t) (T* 


X“ = 0. 


(4.10) 


Here a* are the standard Pauli matrices and hence the matrix size of X* is 4 x 4. Then 


r(t) and x{t) are real fnnctions to be determined. The ansatz (4.10) satishes the Gauss law 
constraint (ig. With the ansatz (|4.10[) , the eqnations of motion (|4.2[) are boiled down to 


0 = r + 2r(r — — 

O 


'■-9 


0 = X + 


2x^ + X + (—2/ir + 6r^) x . 


(4.11) 


One can derive this system (4.11) from the following Lagrangian: 

L = -r-r r-H —x -— x - x — 3r r- ] x 


2 V3 


(4.12) 


Note that the system (4.11) is completely integrable. This means that there exists the 


second integral of motion as well as the Hamiltonian. In fact, the integrability of this system 
was shown in . In order to show the integrability, it is helpful to shift r os r — y -\- jx/Q. 
Then the resulting Lagrangian is 

1 //,a2 


T ^ - 2 . ^ (2 ^ ^ -2 . ^ (2 ^4 q22 

L = -X + X w X — -X + -y + - — y — -y — ox y , 

2 2 V6/ 2 2 2 V6/ ^2 


(4.13) 


where the constant term has been dropped off. Now it is easy to compare the expression 


(|4.13|) with the results in 28 , and the second integral is given by 

2" 


/ = —xy — 2xy 




(4.14) 


It may be interesting to rewrite the Hamiltonian H by subtracting I. One can easily 
show the following expression: 

H-I=^P^-(^^y R^ + ^ Hr . (4.15) 

Here R = x + y and P is the conjugate momentum for R. Thus the Hamiltonian H = Hr + I 
is nothing but a system of two independent oscillators. 
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Finally, it should be remarked that the following orthogonal transformation makes the 
integrability apparent: 


1 

71 



(4.16) 


It recasts X* into the block-diagonal form that indicates that the system is composed of two 
independent pulsating fuzzy spheres with the radial coordinates r-\-x and r — x, respectively. 
Although this ansatz seems rather trivial, it would be worth to mention that the BMN matrix 
model contains an integrable subsystem. 


5 Conclusion and Discussion 


In this paper, we have studied classical chaotic motions in the BMN matrix model. We have 
examined three ansatze: 1) two pulsating fuzzy spheres, 2) a single Coulomb-type potential, 
and 3) two coupled fuzzy spheres. For the first two cases, we have shown the existence of 
chaos by computing Poincare sections and a Lyapunov spectrum. The third case has led to 
an integrable system. As a result, the BMN matrix model is not integrable in the sense of 
Liouville, though there may be some integrable subsectors. 


It would be worth to comment on preceding works on the integrability of the BMN 
matrix model (29] - p^ . The integrability argument is based on the fact that the BMN matrix 
model is obtained as a consistent Kaluza-Klein (KK) reduction of the ^=4 super Yang-Mills 
(SYM) theory on Rx [^. The dilatation operator of the A/'=4 SYM on is identified 
with the Hamiltonian of the J\f=4: SYM on R x via the operator-state mapping. It is 
now realized that the dilatation operator is closely connected with an integrable spin chain 


Hamiltonian [M-36 . Hence the integrability of the dilatation operator might be inherited 
to the BMN matrix model through the KK-reduction. Note here that we have considered 
classical chaos in the BMN matrix model, which is concerned with the classical Hamiltonian 
of the Af=4: SYM on R x (or the classical dilatation operator Dc\ of the Af=4: SYM on 
R^’^). On the other hand, the spin-chain description is realized at the quantum level and 
the quantum dilatation operator plays a central rol^. Thus our analysis has indirect 


^More precisely, is divided as follows: 


-^qu — -^tree 




loop • 


n=l 
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relevance with the preceding one. It may be interesting to study whether the classical chaos 
can survive the quantization. 


The BMN matrix model was originally proposed in the context of M-theory and it can 
be regarded as a deformation of the M(atrix) theory 
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It would be very intriguing to 


try to understand the chaotic motions from the point of view of the DO-brane dynamics. 
With some generalizations, those may capture an aspect of black hole physics such as the 


fast scrambling process 37 . For interesting works along this direction, see 38 . In the very 


recent, a bound on a Lyapunov exponent has been argued 39 


It seems promising that non-linear dynamical perspectives of string theory and M-theory 
would be more and more important. The way and result presented here would be a key 
ingredient in tackling the fundamental problems like information loss due to a black hole 
formation. 
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Appendix 

A How to compute the Lyapunov spectrum 

In Section 2, we have explained the notion of the Lyapunov spectrum in an intuitive way. 
More practically, in order to numerically evaluate the Lyapunov spectrum, we have to follow 
a mathematical formulation. The basic method is based on the expansion rate of the volume 
j^l^. With this method, all of the N Lyapunov exponents can be computed. 

The first term T)tree is concerned with the free Hamiltonian of the BMN matrix model, in which there is no 
chaotic motion as a matter of course. 
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Consider a system of the first differential equations described by 


N\t 


X = F (x) , X = (xy • • • , a:*, • • • , 


(A.l) 


Then, by taking a variation, the following relation is obtained, 

dF^ 

SF = ^SxF 

dx^ 

Hence a deviation vector hx (t) can be expressed as 


(A.2) 


hx (t) = Uf Sx (0) . 


(A.3) 


Here Ut is a time evolution operator and hx (0) is the initial value of the deviation vector. 


This is a solution of (A.2). 


The linearity of Ut leads to the following relation: 


5x {t + s) = Ut5:x. (s) = UfUgS^K (0) . 


(A.4) 


Then the Lyapunov exponent for this deviation vector is defined as 


A = lim log = lim 

T->oo T ||(ix(0)|| n—>-oo nAt 


log 


wuAf.uAtSMm 


(A.5) 


||5x(0)|| 

Here in the last expression T is divided into n small steps like T = nAt . The convergence 
of this limit is proven by Oseledec (4^. 


Let us next introduce a fancy trick to compute the Lyapunov spectrum composed of N 
exponents. The first is to prepare an orthonormal basis set (vj, Vq, ..., v^} to span the 
tangent vector space along the trajectory. Then we make each basis Vq evolve by U/\t ■ The 
resulting set of the vectors {u}, u^, ..., cease to be one of the orthonormal vectors, 
where uj = U^Vq . This is depicted in Fig.[^ 


Then it is necessary to rebuild a new set of the orthogonal vectors {v{, v^,..., from 
{u^, u^,..., by following the Gram-Schmidt process and compute the expansion rate 
for each vector : 


a? = 


= Vi 


Lastly, the orthogonal vector set (v]^, v^,..., is normalized as 

(fc = l,2,...,A). 


Vi = 


1 ~k 

~k^i 

af 


(A.6) 


(A.7) 
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Figure 4: Time evolution of a trajectory and a deviation vector set along it. 


Thus a single cycle has been completed. Then the same cycle should be repeated n times. 
This means that the set {v^^, v|,..., vf^} is used as an initial orthonormal vector set for the 
(z + l)-th cycle. 


Now the largest Lyapunov exponent Ai is defined as 


Ai 


lim 
n-J-oo n/\t 




(A.8) 


i=l 


By construction, vj heads toward the direction which has the highest sensitivity to initial 
conditions at sufficiently late time. Thus the expansion rate aj has the largest value in this 


region. Then it is obvious from (A.5) that the largest Lyapunov exponent is given by the 
sum of logarithm of the expansion rate. 


In a similar way, the k-th Lyapunov exponent is defined as 


Xh = lim 


n 


nAt 


a, 


(A.9) 


1=1 


Thus the set of the exponents {Ai,..., Ayr} is obtained. This set is called the Lyapunov 
spectrum. Note that the computation time can be estimated by the inverse of the smallest 
positive Lyapunov exponent. This is because all of the exponents converge within this scale. 


A characteristic of the Lyapunov spectrum in the Hamiltonian system is that the total 
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sum of the exponents becomes zero, i.e. 


N 




0 . 


(A.IO) 


This relation follows from the computation scheme of the Lyapunov spectrum. By construc¬ 
tion, each of the exponent measures the expansion rate of the orthogonal basis in the phase 
space. Hence the volume V(t) spanned by the bases is expressed as 

V {t) = exp (All) 


Liouville’s theorem says that the phase-space volume is conserved in the Hamiltonian system 


(i.e., V(t) = H(0)). Thus the relation (A.IO) must be satisfied. 

As a side note, for the dissipative system, the following inequality is satisfied: 


N 


Afc < 0. 


(A.12) 


k=l 


This is just because the volume V(t) monotonically decreases. 
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